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j> We study perturbations of exactly tri-bimaximal neutrino mixings under the assumption 
that they are coming solely from the charged lepton mass matrix. This may be plausible in 
\^ scenarios where the mass generation mechanisms of neutrinos and charged leptons/quarks 
^ have a different origin. As a working hypothesis, we assume mass textures which may be 
CN generated by the Froggatt-Nielsen mechanism for the charged lepton and quark sectors, 
which generically leads to strong hierarchies, whereas the neutrino sector is exactly tri- 
bimaximal with a mild (normal) hierarchy. We find that in this approach, deviations from 
^ maximal atmospheric mixing can be introduced without affecting 6'i3 and 612, whereas a 
*^ deviation of ^13 or 612 from its tri-bimaximal value will inevitably lead to a similar-sized 
^ deviation of the other parameter. Therefore, the already very precise knowledge of 612 points 
^ towards small sin^ 2^13 < 0.01. The magnitude of this deviation can be controlled by the 
specific form of the charged lepton texture. 
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1 Introduction 



Comparing the neutrino masses with the charged lepton and quark masses, they observe 
a relatively mild hierarchy. One can easily see that if one expands the masses and mixing 
angles of quarks and charged leptons in terms of powers of a single small expansion parameter 
e. In fact, the CKM mixing matrix Vckm (i||2] exhibits quark mixing angles of the orders 

(1) 



\Vus\ ~ e, \Vcb\ ~ e^ iKfel ~ 



where the quantity e is of the order of the Cabibbo angle e ~ — 0.2. Similarly, for 
the same value e ~ 6'c, the mass ratios of the up-type quarks, down-type quarks, and the 
charged leptons can be approximated, e.g., by 

niu : rric : rrit ^ : : 1, rud : rris : mf, ~ e*^ : : 1, mg : : rrir ~ : : 1, (2) 

where nih/mt ~ e^, m^/mi, ~ 1, and mj ~ 175 GeV. On the other hand, the neutrino mass 
spectrum can be roughly written as 

mi : 7712 : rn^ ~ : e : 1, nii : m2 : ~ 1 : 1 : e, mi : m2 : 7^73 ~ 1 : 1 : 1, (3) 

for the normal hierarchical, inverse hierarchical, and degenerate neutrino mass spectrum, 
respectively. While the mixing angles in the quark sector are small, the lepton Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) mixing matrix t/pMNS (3] exhibits two large mixing angles 
and a small (zero) one ji]. It can be well approximated by the tri-bimaximal (TBM) mixing 
matrix f/xBM fs] (up to phases) as 



PMNS 



TBM 



v 



2 

1 

1 

V6 



1 
1 

1 

V3 



1 

72/ 



35° and 6 



23 



(4) 



45°, whereas 



In f/xBM, the solar and the atmospheric angle are given by 612 
the reactor angle di^ vanishes. 

The experimental (measured) values can arise as deviations from the TBM ansatz of the 
neutrino mass matrix (6]-[9], describing nearly TBM lepton mixing ^lOj. These perturbations 
are often motivated by non-Abelian discrete symmetries 11 , such as A4 
S4 
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The main reason is that there exist suitable breakings of such symmetries into 
subgroups that allow a neutrino mass matrix exactly diagonalized by TBM. To reconcile 
with the experimental data, the perturbations can be induced either by the charged lepton 
sector or by next-to-leading order contributions to the neutrino mass matrix, or both. In 
this work, however, we assume that the mechanism for generating charged fermions and 
neutrino masses and mixings are different. This may not be a remote possibility. In fact, we 
can assume that Majorana neutrino masses are introduced by the lepton number violating 
Weinberg operator [15] , 

Ow = (L^ir^H) (HiT^L) , (5) 

which leads, after electroweak symmetry breaking (EWSB), to Majorana masses for the 
neutrinos {H is the SM Higgs doublet). It is well known that this operator implies physics 
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beyond the Standard Model, such as a heavy neutral fermion leading to the type I see-saw. 
Therefore, the origin of the Majorana neutrino mass lies in physics beyond the Standard 
Model, including couplings beyond the Standard Model. On the other hand, the other 
fermion masses can be easily described within the Standard Model, at least as long as the 
hierarchies need not to be justified. Or, invoking an 5'[/(5)-inspired grand unified framework, 
they can be understood because charged leptons and quarks are arranged in the same GUT 
multiplets (the neutrinos being singlet of the group). 

Therefore, it may be natural to assume that the leading mass generation mechanisms of 
neutrinos versus charged leptons/quarks are different. There is, however, one drawback 
of this strictly phenomenological separation: since neutrinos and charged leptons come in 
SU(2) doublets in the Standard Model, this strict separation above the EWSB scale might 
be challenging. For example, left-handed neutrinos and charged leptons will always belong 
to the same representation of the non-Abelian discrete group. On the other hand, one or 
both of the mass generation mechanisms of the neutrinos and charged leptons/quarks may 
be implemented at or below the EWSB scale. We will not enter this level of detail, but 
instead study the phenomenological consequences of such a model with mass generation 
mechanisms, which are separated to leading order, can be constructed. 

As one example for the charged lepton and quark mass generation, we use the Froggatt- 



Nielsen (FN) mechanism 16 in which effective dimension-n mass terms lead to masses 
proportional to e", where e depends on the flavon vacuum expectation value suppressed 
by the mass of super-heavy fermions. In this way, mass matrix textures with e-powers as 
entries are obtained. Consequently, such a matrix structure contains information on the 
hierarchy among matrix elements and goes beyond approaches which use texture zeros. 
The FN mechanism is a perfectly plausible possibility to generate strong hierarchies. It has, 
for instance, been used in Refs. Il7j to construct charged lepton and even neutrino mass 



textures which can be implemented by discrete flavor symmetries 13 . Here we use a similar 
approach to study the interplay between the charged lepton mass matrix, generated by this 
approach, with TBM mixings in the neutrino sector and small deviations from it coming 
from the diagonalization of the charged lepton mass matrix. Note that in this ansatz the 
quantity e determines both the charged lepton and quark mass matrices, which points to a 
common origin, leading to some form of "quark-lepton complementarity". QLC has been 
studied from many different points of view 18 -24 . As in the earlier references 17 , we 



implement this quark-lepton complementarity at the Yukawa coupling level, but we study 
the implications of random complex order one coefficients, as suggested by the original 
Froggatt-Nielsen approach, within specific textures leading to interesting deviations from 
TBM mixings. 

2 Methods 

We diagonalize the charged lepton and Majorana neutrino mass matrices as 

Me = U,mY'"^U[^ (6) 
= U^^C^^U^, (7) 
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Best-fit 


3(7 (current) A 


3a (2015) A 


3a (2025) A 


3a (2035) A 


sin^^i2 0.318 
sin^^ig 0.013 
sin2^23 0.5 


0.27. ..0.38 

< 0.053 e 
0.36. ..0.67 e 


No furt 
< 0.012 e 
0.43. . .0.57 


her experiments f 

< 0.001 
0.47. . .0.53 


)lanned? 
< 1.5 ■ 10-5 
0.47... 0.53 



Table 1: Current best-fit values for the mixings angles and the 3cr allowed ranges '4h well as projections 



for the measurements labeled "2015" (mostly Daya Bay and T2K, see Ref. 120^ ), "2025" (superbeam up- 
grades, such as T2HK or LBNE; see Ref ^ for sin^ 6*13 and Ref ^ for sin^ 623), and "2035" (neutrino 
factory '29,30^); for the projections, sin^ ^13 = is assumed as best-fit value. The columns labeled "A" 



motivate typical allowed deviations of the angles from their TBM values in terms of the e power such that 
A is comfortably with the 3(t allowed range for order one coefficients A in Eq. (10) (for e = 0.2). 



where M°'^^ and M°'^g are, up to an overall mass scale, given by Eq. ^ (third rela- 
tionship) and Eq. (|3|, respectively. The matrices Ui, U[, and U^, are in general, arbitrary 
unitary matrices. By our assumptions, we choose = U^bm = U{sm'^ 612 = 1/3, sin^ ^13 = 



25 



There- 



0,sin^6'23 = 1/2) with U the mixing matrix in the standard parameterization 
fore, the neutrino mass matrix is, together with the choice of the hierarchy in Eq. ^ and 
the absolute neutrino mass scale, uniquely determined and assumes special versions of the 
TBM form. The lepton mixing matrix is, as usual, given by 



PMNS 



(8) 



Using Eq. ^ in Eq. ([6]), we obtain 

Me = ([/TBMt/fMNs)^?" 'f^? > (9) 

where ?7pmns is the measured mixing matrix. Note that Eq. ^ allows us to construct Mi 
for arbitrary t/pMNS deviating from TBM mixings if we fix U^. In the following, we will use 
U[ = 1 for the sake of simplicity, which means that = Mf^'^^ if t/pMNS has exactly the 
TBM form. On the other hand, any deviation from TBM mixings will lead to a non-trivial 
mass matrix for Mi with off-diagonal entries. Our approach therefore corresponds to a 
particular form of perturbations of the TBM mixings coming from the charged lepton mass 
matrix only. 

Since we specify the mass hierarchies in terms of e, we will also motivate the deviations 
from TBM mixings by powers of e. We define 



Oi'i — 9, 



TBM 



{12,13,23} 



(10) 



with order one coefficients e < \Aij\ < 1/e. We list in Table [T] the estimated experimental 
ranges for four different experiment generations, denoted by "Current" (current best-fit), 
"2015" (mostly from Daya Bay and T2K), "2025" (superbeam upgrades), and "2035" (neu- 
trino factory). In this table, we also list plausible deviations in terms of the power of e from 
TBM mixing for each angle, motivated by these measurement precisions. One can clearly 
read off the table, that each generation of experiments will improve the precision on the 
mixing angles. If the TBM values are confirmed, Ajj will be constrained further, and smaller 
deviations from TBM mixings will be allowed. Very importantly, sin^ 612 is already very well 
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measured oc e^, which has interesting imphcations - as we will demonstrate later. We will 

as input assumptions instead of the measurement precisions, 



further on use A 



but Table [T] demonstrates that these are closely related. Using this hypothesis for t/pMNS in 
Eq. ([9]), we can construct Mi for each case. In that case, the mass matrix Mi depends on e 
(and the absolute mass scale) only. In the spirit of Ref. 17 , we can then extract a texture 



from the mass matrix by identifying the leading entry and absorbing the lowest power of e 
in the absolute neutrino mass scale. For example, we identify (mass matrix — )■ texture): 



:iii 



where "0" stands for C(e^). 

In the reverse direction, specifying as given by the theory and re-diagonalizing it, we 
should be able, barring ambiguities, to reconstruct the initial hypothesis for the deviation 
from TBM mixingsj^ For example, the texture can then be interpreted in terms of Froggatt- 
Nielsen-like models with arbitrary (complex) order one coefficients Cij: 
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:i2) 



In the Froggatt-Nielsen ansatz, the order of the texture entries is given by a discrete symme- 
try, whereas the order one coefficients Cjj are random order one entries. Where applicable to 
test the stability of individual charged lepton textures, we generate these order one entries 
Cij = \cij\ exp{i(j)ij) randomly, with |cjj| uniformly distributed between e and 1/e and (pij in 
[0,27r[F] Once M^, given by Eq. (g, and M^, given by Eq. ([t]), are specified, the two mass 
matrices can be numerically diagonalized, ?7pmns can be computed, and the mixing param- 
eters can be read off from C/pmns using re-phasing invariants. We use the MPT package 



from Ref. 31 for this part. 



3 Deviations coming from charged lepton mass matrix 



First, we investigate the effect of arbitrary (small) deviations from TBM mixing coming from 
the charged lepton mass matrix. Our notation (Aj^) can be easily related to other notations 
used in the literature. Here we use the parameterization of deviations from TBM from 
Ref. p], which is in terms of the sines of the angles, not the angles themselves. Therefore, 
additional pre-factors appear: 

a ~ A23 , s ~ V2A12 , r ^ V2A13 . (13) 

^Note that U' is, by definition, not relevant for this method, because the re-diagonahzation will yield 
U' ~ 1 automatically if the mass matrix is constructed that way. 

^In fact, we draw k between e and 1 uniformly, and then assign \cij\ — k ot \cij\ — 1/k with 50% 
probability each. 
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We use Eq. (13) in Eq. ^ and expand the resulting charged lepton mass matrix to second 
order. In this way, we obtain 



M: 



general 



/ 1 



v 



r+s 

2 
r~s 

2 



1 - 

a 



r{l+a)+s{l-a) 

_ (r+s)^/2+2a^ 
4 



1 



s—r+{r+s)a 
2 

_ (r+s)^/2-s^ 

' 4 

_ (r-s)^/2+2a^ 

4 



diag(me,m^,mr) 



(14) 



This means that our approach lead to the universal mass matrix Eq. (14) which is solely 
determined by observables. Moreover, the one-to-one correspondence of Mf^"^^"^^^ with the 



experimentally accessible quantities in Eq. (13) in fact allows us to determine the Yukawa 



couplings in this approach, which can be extracted using Eq. ( 14 ) 



From the general form of Mf'^'^'^'^^^ in Eq. (14), we can already make some interesting obser- 



vations. Suppose we wanted to have individual deviations from TBM mixings, such as we 
only wanted deviations in the 2-3 sector, i.e., r = s = 0: 



1 
1 - 








a 



91. 
2 



diag(me, m^,m^) 



(15) 



Using Eq. (|8]), it is easy to verify that the first row of ?7pmns, which determines sin6'i3 and 
tan6'i2, remains fixed to the corresponding TBM entries. This means that we can introduce 
deviations to ^23 = 7r/4 without affecting ^13 and 9i2 in this approach. In the reverse 
direction, the conclusion is that measuring ^13 = together with a deviation from maximal 
atmospheric mixing can be trivially interpreted as a perturbation from the charged lepton 
sector. 



On the other hand, if we choose a = 0, we obtain from Eq. (14): 



m; 



c=0 



/ 1 

v 



r+s 

2 
r—s 

2 



1 - 



r+s 
(r+sf 



-2rs+s^ 



—r'^+2rs+s^ 



Ir-sf 



\ 



J 



diag(me, m^,m^) 



(16) 



Here we see that the matrix is symmetric in deviation of ^^13 and 9i2- The consequence will 
be that a perturbation in one of the sectors affects also the other mixing angle, because if it 
is introduced at the mass matrix level, it cannot be clearly assigned to one of the angles. For 
example, this can also be seen very specifically in the special case of deviations introduced 
in the 12-sector of the charged lepton mixing only (CKM-like charged lepton mixings). The 
sum rule from Refs. 32 reads, in our notation, A12 — A13 cos 5, which shows that the two 



mixing angles are related. 



4 Mass textures for deviations from TBM mixings 

In this section, we discuss how we can derive charged lepton mass textures in terms of 
powers of e for certain perturbations of TBM. We will discuss in the next section, if these 
textures, if given by a model, indeed lead to the expected deviations from TBM using 
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the Froggatt-Nielsen ansatz as example. We here parameterize the deviations from TBM 
mixings in terms of e: 



A 



13 



Re 



A 



12 



"12 



A 



23 



Ae 



"23 



(17) 



where riij 



1, 2, 3, . . . is the power of e which determines the (leading) magnitude of the 
deviation, and R, S, A are order one coefficients, i.e., e < \R\, \S\, \A\ < 1/e. We can also 
parameterize the charged lepton and neutrino masses in terms of e, see Eg. pi and Eq. (§, 
where we use the normal hierarchy case as an example. Applying Eq. (17) and the mass 



spectra to Eq. (14), we can derive mass textures for individual perturbations of TBM in 



terms of e-powers; c/., Eq. (11) 



As a first example, let us introduce a deviation from maximal atmospheric mixing assuming 
R = S = 0. From Eq. (15), we read off that for A23 = Ae^ 



c"23 



Mf 








^ 


















"23 = 1 j 












^g(2+n23) 


1 1 




^ 


e3 





:i8) 



where we have computed the texture for n23 = 1 in the last step. As it is obvious from the 
discussion in the previous section, the choice of this texture does not significantly affect the 
TBM values of ^13 and 9i2. We have checked that this conclusion even holds in the presence 



of arbitrary order one coefficients in the individual texture entries; c/., Eq. (12). The size 



of the deviation from TBM is given by the texture entries in Eq. (18), i.e., by choosing a 
particular texture, one can control the magnitude of the deviation from TBM. 

If we want to introduce deviations to the TBM values of ^13 or 9i2, we know from Eq. (16) 



that these can be not easily generated separately. In general, we find that the texture 
is determined by the leading (largest) deviation from TBM, or the lowest power n = 
min(ni3,ni2) > 1, as 








g2+n 
,2+2n 



^2n 



(19) 



This leads for n 



1, 2, or 3, to three distinct textures: 









e 













V 




1 



n=2 













(20) 



It is characteristic for that texture that the smaller the deviation from TBM is, the more 
so-called texture zeros [(9(e^)] are in the charged lepton textures. Moreover, the structure 
as a whole, represented by powers of e in the matrix elements, does not change, i.e., the e 
power of each matrix element ascends or remains constant while going from case n = 1 to 
n = 3. For the particular case n = ni2 = rii^, the mass matrix reads explicitly to leading 
order in e 



Mf 





V 



A. 

'V2 



Ah 



2+n 



-R^ + 2RS + S^) e2+2- 



( _S R_\ 

_ 2RS + S^) e^" 
1 



(21) 
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The case ni2 > nis can be obtained easily from Eq. (21 ) by setting S = 0, the case > ni2 
by setting R = 0. Here the 0{1) coefficients are exphcitly taken into account through the R 
parameter. This mass matrix allows us to interpret future possible measurements of di^ and 
6i2 in terms of the three different models ?T, = l,n = 2, orr;, = 3, which may be generated 
by discrete flavor symmetries, i.e., the Yukawa couplings can be measured. 

5 Stability of approach in Froggatt-Nielsen models 



In Froggatt-Nielsen models 16 , the Yukawa couplings may arise from higher-dimension 



terms in combination with a flavor symmetry: 

Ces-^{H)e^^L^R. (22) 

In this case, e becomes meaningful in terms of a small parameter e = v/Mp which controls 
the flavor symmetry breaking]^ The integer power of e is solely determined by the quantum 
numbers of the fermions under the flavor symmetry. Therefore, the flavor symmetry pre- 
dicts the mass texture at the Yukawa coupling level, where, in the original FN approach, 
the coefficients are arbitrarily chosen complex order one numbers. We therefore test the 
stability of our textures in this framework using random order one coefficients, as described 



in Eq. (12) 



Introducing deviations from 6*23 = 7r/4 only, i.e., by Eq. (18), does not produce any qualita- 
tively new insight: using this texture, the other two mixing angles are basically not affected. 
This means that a deviation A23 can indeed be introduced in this model without perturbing 
the TBM values of ^13 and 612- 



More interesting are the textures in Eq. (|20j). We show the effects from constrained variation 
of Yukawa couplings for n = 1 (left column), n = 2 (middle column), and n = 3 (right 
column) in Fig. [l] Here we can immediately see the implication of Eq. ( [2l| ): since the leading 
order entry in each mass matrix element dominates, the order one coefficients will induce 



both Ai3 and A12 with the same (leading) order in Eq. (17), i.e., ni3 = ni2. Therefore, 
for arbitrarily chosen 0{1) coefficients, a large 613 comes together with a large deviation of 
9i2 from its TBM value. In addition, similar-sized deviations of 623 from maximality are 
introduced by these charged lepton textures. 

In Fig. [1] we show the current bounds and impact of future experiment generations if the 
TBM values are conffimed. While the precision of ^13 and 623 do currently not put the 
case n = 1 under pressure, the high precision to which 612 is measured only allows few 
possibilities for the n = 1 case. This is illustrated by the first column of Table [l| where it 
is indicated that A12 oc e^, whereas A13 and A23 oc e. Therefore, A12 is the discriminator 
in this case, and n = 1 can be basically excluded, and with it, large values of ^13, i.e., 
sin^ 2^13 < 0.01, as we can read off from the upper middle panel. This situation changes for 
future experiments, where the more precise measurements of ^13 and 623 will constrain the 
textures further: the cases "2015" and "2025" enter the parameter space of n = 2, whereas 



^ Here v are universal VEVs of SM singlet scalar "flavons" that break the flavor symmetry, and Mp 
refers to the mass of super-heavy fermions, which are charged under the flavor symmetry. The SM fermions 
are given by the ^'s. 
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Figure 1: Correlations among the PMNS mixing angles for 600 unconstrained variation 
of Yukawa couplings of Eq. (20) for n = 1 (left column), n = 2 (middle column), and 
n = 3 (right column). In any panel, the Yukawas are complex parameters with moduli 
between 1/e and e and unconstrained phases in [0,2tt[, where e = 0.2. The shaded regions 
are experimentally excluded and the sensitivity reach of future neutrino experiments on the 
various mixing angles are marked (the current 3a bounds, the expected limit in "2015" 
(dashed line), the limit in "2025" (dash-dotted line) and the reach in "2035" (dotted line)). 
The circle indicates the best fit points for the mixing angles, according to Table [I[ 
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the Neutrino Factory may even constrain n = 3 by its extremely well sin^ 26'i3 bound. As one 
can read off from the right panels, the pressure on it will then indirectly also constrain 
the deviations of 9i2 and 6*23 from their TBM values. 

6 Summary and conclusions 

We have illustrated that introducing deviations to TBM from the charged lepton sector 
only is a viable and simple approach, which may be motivated by models where the mass 
generation mechanisms of neutrinos and charged leptons/quarks decouple. Deviations of 623 
can be independently induced by a particular charged lepton texture without affecting the 
TBM values of ^13 and 612- The specific form of the texture then controls the magnitude 
of the deviation from TBM. In a similar way, deviations of ^13 and 612 can be introduced, 
where the specific form of the texture again controls the magnitude of the deviations. In 
this case, however, the deviations of all parameters from TBM are typically introduced with 
the same magnitude. We have tested and confirmed these conclusions in a Froggatt-Nielsen 
approach for the charged lepton mass matrix, where we have chosen order one (random) 
coefficients. In addition, we have shown that the deviations from TBM can be used in 
particular models to directly measure the Yukawa couplings. 

As the main conclusion in this approach, the entanglement between ^13 and 612 implies that 
the currently extremely good measurement of 6*12 already exerts pressure on 613. In fact, 
plausible textures are obtained for sin^ 2^13 < 0.01. On the other hand, future measurements 
of di3 will limit the deviations of 612 and 6*23 from their TBM values. Deviations of ^23 = 7r/4 
can, however, be independently introduced, which means that any measured such deviation 
will not lead to any conclusions for 613 and 612 in this approach. 
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